arXiv: 1504.07712vl [cond-mat.stat-mech] 29 Apr 2015 


Critical behaviors and universality classes of percolation phase transitions on 

two-dimensional square lattice 

Yong Zhu/ Ziqing Yang/ Xin Zhang/ and Xiaosong Chen^ 

^ State Key Laboratory of Theoretical Physics, Institute of Theoretical Physics, 

Chinese Academy of Sciences, P.O. Box 2735, Beijing 100190, China 
(Dated: April 30, 2015) 

We have investigated both site and bond percolation on two dimensional lattice under the random 
rule and the product rule respectively. With the random rule, sites or bonds are added randomly 
into the lattice. From two candidates picked randomly, the site or bond with the smaller size product 
of two connected clusters is added when the product rule is taken. Not only the size of the largest 
cluster but also its size jump are studied to characterize the universality class of percolation. The 
finite-size scaling forms of giant cluster size and size jump are proposed and used to determine the 
critical exponents of percolation from Monte Carlo data. It is found that the critical exponents of 
both size and size jump in random site percolation are equal to that in random bond percolation. 

With the random rule, site and bond percolation belong to the same universality class. We obtain 
the critical exponents of the site percolation under the product rule, which are different from that 
of both random percolation and the bond percolation under the product rule. The universality class 
of site percolation differs different from that of bond percolation when the product rule is used. 
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I. INTRODUCTION 

Traditionally, percolation phase transitions are con¬ 
sidered to be continuous in various networks. How¬ 
ever, Achlioptas et al. [I| concluded that the percola¬ 
tion transition in random network is discontinuous un¬ 
der the Achlioptas process (AP) in which the growth of 
large clusters is suppressed. During the evolution of a 
network under AP, edges are added into the network un¬ 
der product rule that the edge with minimum product of 
cluster sizes is connected from two randomly picked unoc¬ 
cupied edges. The percolation transitions under product 
rule in two-dimensional square lattice and scale- 

free networks [3-Q are also investigated and argued to 
be discontinuous. Later, a series of investigations 0- 
[TTI . [T^ showed that different percolation transitions un¬ 
der product rule are actually continuous. In Ref. Q, 
Costa et al. Q studied a model with stronger suppression 
of large clusters and its percolation transition is contin¬ 
uous. Lee et al. Q concluded that the explosive perco¬ 
lation is continuous by studying cluster size distribution. 
Grassberger et al. @ found that the explosive transition 
is still continuous but with unusual finite-size behavior. 
In the works of Liu et al. [l^ and Fan et al. , a gen¬ 
eralized Achlioptas process (GAP) was introduced to in¬ 
vestigate the percolation transitions of two-dimensional 
lattice network and Erdos-Renyi (ER) network [T^ re¬ 
spectively. In the generalized Achlioptas process, the 
edge with minimum product of cluster sizes is connected 
with a probability p from two randomly picked unoccu¬ 
pied edges. When p = 1/2, GAP recovers to a random 
process. With p = 1, GAP becomes AP. The percola¬ 
tion transitions of two-dimensional lattice network and 
ER network under GAP keep to be continuous. Their 
critical exponents and and therefore universality classes 


depend on the probability parameter p. With a rigorous 
mathematical proof, Riordan and Warnke [l^ concluded 
that the percolation transitions under all Achlioptas pro¬ 
cesses are continuous. 

The site percolation on a two-dimensional square lat¬ 
tice under product rule was first investigated by Ghoi et 
al. H- They claimed that the percolation is discontin¬ 
uous from the non-vanishing hysteresis between the di¬ 
rected and reverse process for the fraction of the sites be¬ 
longing to the largest cluster. However, Bastas et al. [l^ 
demonstrated that the hysteresis phenomena will disap¬ 
pear in the thermodynamic limit. Using the finite-size 
behaviors of the largest cluster and its standard devia¬ 
tion, they concluded that the explosive percolation tran¬ 
sition is continuous but belongs to different universality 
class. To clarify this controversy, further investigations 
are needed. 

The notation of a universality class is a basic tenet in 
the physics of critical phenomena. Within a universality 
class, the universal quantities such as critical exponents 
and scaling functions are independent of microscopic de¬ 
tails. It is accepted traditionally that a universality class 
is characterized by the spatial dimensionality of the sys¬ 
tem and by the number of the components of the order 
parameter. (See, e.g., the review article [l^.) On a two- 
dimensional lattice, the random bond percolation has the 
same critical exponents as the random site percolation. 
The random bond and site percolation belong to the same 
universality class. For the bond percolation under GAP 
with p > 1/2, their critical exponents are different from 
that of random bond percolation and depend on the prob¬ 
ability parameter p [l^. The universality classes of bond 
percolation under GAP are different from that of random 
bond percolation 0- It is of interest to investigate the 
critical behaviors and universality class of site percola- 


2 


tion in relation to the Achlioptas process. 

In this paper, we study the critical behaviors of site 
percolation under AP in two-dimensional square lattice. 
We investigate the sizes of giant clusters in the lattice. 
From the finite-size scaling behaviors of giant clusters, 
we can determine the percolation transition point and 
the corresponding critical exponents. Both the critical 
exponents of the sizes of giant clusters and their size gaps 
during evolution [^, [13, E] are calculated. 

Our paper is organized as following. In Section 2, we 
present the finite-size scaling behaviors of giant clusters 
near the phase transition point of percolation. The finite- 
size scaling behaviors of giant clusters for site percolation 
under random and product rules on two-dimensional lat¬ 
tice are studied in Sections 3 and 4 respectively. The 
universality classes of bond and site percolation under 
random and product rules are discussed in Section 5. 


II. CRITICAL BEHAVIORS OF PERCOLATION 
PHASE TRANSITION 


which is a straight line against InL. We can fix the criti¬ 
cal point Tc by the straight line of In si (r, L) with respect 
to InL. From the slope of the straight line, we can de¬ 
termine the critical exponent ratio fijv. By introducing 
scaling variable to make curves s\ of different 

L collapse into a scaling function the critical 

exponent u can be estimated. 

There are critical behaviors also during the evolution 
of a network. The percolation phase transition is accom¬ 
panied not only by the appearance of a giant cluster but 
also by the giant size jump of the largest cluster. In the 
Refs. 171 and [l^, the percolation phase transition was 
investigated by the size jumps of the largest cluster. 

In a simulation of network evolution, edges are added 
under some rule one by one into network. At an evo¬ 
lution step T, the number of edges in the network in¬ 
creases from T — 1 to T and the size of the largest clus¬ 
ter varies correspondingly from S'i(r — 1) to S'i('L). At 
this step, the largest cluster has a reduced size jump 
St = [«S'i(L) — Si(T — 1)] /N. The largest reduced size 
jump of the whole evolution process is 


The percolation phase transition in a network with 
N = L‘^ nodes is indicated by the appearance of a giant 
cluster whose size becomes comparable with N. The size 
Si of the giant cluster is taken as the order-parameter of 
percolation phase transition. Near a critical point, the 
size S 2 of the second largest cluster in the network be¬ 
comes comparable with N also. Both Si and S 2 demon¬ 
strate critical behavior near the critical point llll|. 

Following the finite-size scaling of order-parameter |23l . 
1^ . we anticipate that the reduced sizes of and S 2 in a 
network with N nodes and Nr edges follow the finite-size 
scaling form as Elm 


siir,L)=Si/L‘^ = SiitL^^n. ( 1 ) 

S2{r,L)=S2/L<^ = L-P/''S2{tL^'''). ( 2 ) 


where r = Nr/N is the reduced number of edges and 
t = — is the deviation from critical point Va- The 

scaling variable is related to the size ratio L/^ of L 
to the correlation length in bulk ^ The finite- 

size scaling form is supposed to be valid in the asymptotic 
critical region with L ^ 1 and |t| ^ 1. 

From Eqs. © and @, the finite-size scaling form of 
size ratio S 2 /S 1 is obtained as 

S 2 /S 1 = = U{tL^/n- (3) 


At the critical point Tc, the size ratio S 2 /S 1 = U{0) 
and becomes independent of system size L. The curves 
S 2 /S 1 of different L against r have a cross-point, which 
can be used to determine the critical point re- 
The logarithm of Eq. is 

lnsi(r, L) = —(P/v) InL -|- lnsi(fL^/‘^). (4) 


(5) 


A = max {<51,^2, (6) 


The evolution step, where the largest size jump ap¬ 
pears, is denoted as T^. In the Lth simulation, we obtain 
the largest reduced size jump and the critical re¬ 
duced evolution step Vc'^ = Tc'’ /N. Erom the results of 
M simulations, we can calculate following averages 


A(l) 


rc(L) 


1 

M 

1 

M 


M 

i^l 

M 

(7) 

2=1 

(8) 


In the bulk limit L ^ 00 , we suppose that rc{L) ap¬ 
proaches its bulk limit rdoo) as 


rc(L) = rdoo) + UrL (9) 


The bulk limit re ( 00 ) from the network evolution should 
be equal to the critical point Tc obtained from the largest 
cluster. 

The character of percolation phase transition can be 
determined from the finite-size effect of A(L). If A(L) 
approaches a non-zero value in the limit L —>■ 00 , the per¬ 
colation is a discontinuous phase transition. Eor a contin¬ 
uous percolation phase transition, A(L) has a power-law 
finite-size effect as 


A{L)=aAL-^\ (10) 

In the i-th simulation of network evolution, there are 
fluctuations Sre = r^'’ — fc{L) and dA = A^®^ — A(L). 
Their root mean squares are calculated as 

Xr = \/< {Sre)"^ >, 

XA ^ V< (<5A)2 >. 


At the critical point, we have 

lnsi(rc, L) = —(/3/i/) InL -f Insi(O) 


( 11 ) 

( 12 ) 
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Their dependence on network size L is described by 
the exponents i >2 and /32 as 

Xr OC (13) 

XA OC L~^^. (14) 

In Ref. [13 , the asymptotic behavior of A(L) was used 
to judge the continuity of a percolation phase transition. 
The so-called upper pseudo transition point in Ref. @ 
was estimated from rdL)- The general clique percolation 
phase transition in random networks are identified by 
studying the averages and fluctuations of A and Tc [T^ - 
We will investigate the above critical behaviors of a 
percolation phase transition, which are characterized by 
the critical exponents /3, Pi,132, v,vi and V 2 . The differ¬ 
ent percolation phase transitions with the same critical 
exponents belong to the same universality class. 


III. RANDOM SITE PERCOLATION ON 
TWO-DIMENSIONAL SQUARE LATTICE 

On a two-dimensional square lattice, the percolation 
transition can be investigated both for site and bond. For 
the site percolation, we consider a L x L square lattices 
with periodic boundary conditions in both directions. At 
first, there are N = L x L unoccupied sites in a lattice. 
Then the sites in the lattice are occupied one by one 
under some rule. Clusters consisting of adjacent sites 
appear during the evolution. When Nr sites are occupied 
in the lattice, we define a reduced number of occupied 
sites r = Nr/N. If the size of the largest cluster becomes 
comparable with N, there is a site percolation in the 
lattice. 

For the random site percolation (SP), sites are added 
randomly. It is known that SP on two-dimensional 
square lattice occurs at Tc = 0.59276421(13) [l^. From 
Coulomb gas arguments and conformal field the¬ 
ory [2l|, the critical exponents P = 5/36 and = 4/3 
were predicted. 

In our Monte Carlo simulation of site percola¬ 
tion, the algorithm of Newmann and Ziff |T^ . 
has been used. We have taken lattice sizes L = 
32, 64,128, 256, 512,1024,2048,4096 and 8192 in our sim¬ 
ulations. 12, 800,000 independent simulations have been 
run for each lattice size. 

In FiglU the size ratio S 2 /S 1 is plotted with respect to 
r for three different lattice sizes. From its cross-point, we 
get the critical point Tc = 0.59276(3) which agrees with 
the results of Ref. [l^ . 

In Figl^ In Si is shown with respect to InL for three 
different r. At r = 0.5927, the curve is nearly a straight 
line with slope equal to —0.108. The curvature becomes 
negative for r = 0.5924 and positive for r = 0.5930. 
We can conclude that the critical point Tc = 0.5927(3) 
and the critical exponent ratio Pjv = 0.108(4), which 
agrees well with the exact value Pjv = 5/48. Within er¬ 
ror bounds, the critical point from In si agrees with that 



Figure 1: Size ratio of the second largest cluster to the largest 
cluster U = S 2 /S 1 . The critical point of site percolation = 
0.59276(3) is obtained from its cross-point. 



Figure 2: Log-log plot of si with respect to L for three differ¬ 
ent reduced number of occupied sites r. At = 0.5927(3), 
the curve becomes a straight line and its slope gives the crit¬ 
ical exponent ratio Pjv = 0.108. 


from S 2 /S 1 . Using Ijv = 0.75 in FiglSj three curves of 
S 2 /Si at L = 256, 512,1024 collapse into a scaling func¬ 
tion function 

In Figdl we show critical behaviors of the largest clus¬ 
ter during the evolution of SP. The Monte Carlo sim¬ 
ulation results of fc{L) are demonstrated in Fig|4l(A). 
Using rc(oo) = 0.5927(3) from the above investigation, a 
log-log plot of rc(oo) — fc{L) with respect to L is made. 
Our Monte Carlo data demonstrate a linear dependence 
of ln[rc(oo) — fc{L)] on InL. From its slope, we obtain 
the inverse of exponent tjvi = 0.75(1). In Fig0KB), 
the power-law finite-size behavior of A(L) in Ea. dTUl) is 



,l1/v 


Figure 3: Finite-size scaling function U(tL^^‘') of the size ratio 
S 2 /S 1 for random site percolation with Xjv = 0.75. 
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ln(L) 


ln(L) 


Figure 4: Critical behaviors of random site percolation dur¬ 
ing network evolution: (A)log-log plot of rc(oo) — fc{L) with 
rc(oo) = 0.5927 and the exponent Xjvx = 0.75(1); (B)log-log 
plot of A(Z/) with the exponent /32 = 0.104(1); (C)log-log plot 
of Xr with the exponent 1/(^2 = 0.74(1); (D)log-log plot of xa 
with the exponent — 0.104(1). 


confirmed by the Monte Carlo simulation data. The ex¬ 
ponent /3i = 0.104(1) is calculated from the slope of the 
straight line. 

The root mean squares of fluctuations Xr and xa are 
shown in FigsHKC) and (D). Our Monte Carlo simulation 
data of Xr follow a power-law dependence on L with the 
exponent I/ 1/2 = 0.74(1). The power-law behavior of xa 
in Eg 1 141 is confirm by our simulation data and we obtain 
h = 0.104(1). 

We have investigated also the evolution critical be¬ 
haviors of random bond percolation (BP) on two- 
dimensional square lattice. These results of BP are given 
in TablelH We can see that the critical exponents of SP 
are equal to that of BP. Further, the critical exponent 
ratio fijv is equal to /3i and P 2 and u is equal to ui and 
1^2 within error bounds. 


IV. SITE PERCOLATION UNDER PRODUCT 
RULE ON TWO-DIMENSIONAL SQUARE 
LATTICE 



r 


Figure 5: Reduced size of the largest cluster as a function of 
r at lattice size L = 512, 1024, 2048, 4096, 8192. 



Figure 6: Size ratio of the second largest cluster to the largest 
cluster U = S 2 /S 1 at different lattice sizes. Its fixed point 
gives the critical point = 0.7634(2). 


ical point Tc = 0.7634(5), which is in agreement with 
the value of the fixed point. The accuracy of our critical 
point is higher than that of Tc = 0.768(3) in Ref. 0 and 
rc = 0.756(6) in Ref.0. 

The slope of the straight line in Figl3 gives the critical 
exponent ratio P/u = 0.010(1), which is different from 
Pjv = 0.04 ± 0.02 of Ref.0. 

We can see in Figl8]that the size ratios S 2 /S 1 of dif- 


In this section, we consider the site percolation under 
product rule (SPPR). We take the product rule 00 
for the site evolution as follows: (l)selecting two unoccu¬ 
pied sites randomly; (2)calculating the product of sizes 
of the clusters which are connected by two chosen sites 
respectively; (3)the site with smaller product is occupied. 

The reduced sizes si of the largest cluster for lattice 
size L from 512 to 8192 are shown in FiglSJ si is nearly 
zero at small r and becomes finite when r is larger than 
a critical value r^. Finite si indicates the emergence of a 
giant cluster and a site percolation transition. 

The results obtained from the size of the largest clus¬ 
ter are shown in Figs. [SI 13 and [51 From the fixed 
point of S 2 /S 1 in FigjSl we can get the critical point 
Tc = 0.7634(2). The straight line in Fig|3 gives the crit- 



ln(L) 

Figure 7: Insi with respect to InL around the critical point 
Cc = 0.7634(5), where the curve becomes a straight line with 
slope equal to —0.010. 
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Figure 8: Finite-size scaling function U{tL^^'') of size ratio 
S 2 /S 1 with \/v = 0.42. 


ferent lattice sizes collapse well with the scaling variable 
tV-!'' for data of L > 512. This means that the asymp¬ 
totic region of SPPR is more restricted than that of SP. 
In the finite-size scaling plot of S 2 /S 1 , we use the inverse 
of the critical exponent Xjv = 0.42(1). 

The size ratio at critical point 1/(0) is supposed to be 
universal. We obtain t/(0) = 0.03 for SPPR and t/(0) = 
0.1 for SP. This difference of C/(0) indicates also that 
SPPR and SP belong to different universality class. 

We have studied also the size jumps of the largest 
cluster during network evolution for SPPR. The differ¬ 
ent quantities related to the size jump are shown in 
FigUlwith respect to the lattice size L. From Figl^A), 
we can conclude that the average transition point ap¬ 
proaches to its bulk limit according the power-law in 
EgfTTl Using the Monte Carlo data of the lattice size 
I > 512, we get the inverse of exponent Xjvi — 0.42(2). 
The root mean square of the fluctuations of transition 
point is plotted in FiglOKC) and a power-law in EafT^is 
found l/zz 2 = 0.42(3). The simulation results of the size 
jumps of the largest cluster are presented in FigsEKB) 
and (D). The average size jump A shows a quite differ¬ 
ent L-dependence from that of the random site percola¬ 
tion. No power-law behavior is found for A(L), which 
increases with L at first and then decreases. The root 
mean square of fluctuations of A follows a power-law of 
the lattice size also. From Monte Carlo simulation data, 
we get P 2 = 0.008(4). 

For comparison, we have studied also the evolution 
critical behaviors of bond percolation under product rule 
(BPPR). The critical exponents related to the cluster 
size jump /3i = 0.044(5) and (32 = 0.041(4) are ob¬ 
tained. For the exponents related to critical point, we 
get l/vi = 0.99(5) and l/i /2 = 0.94(3). The critical expo¬ 
nent ratio fi jv and the inverse of critical exponent 1/zz of 
BPPR have been calculated in Ref.[l3. fi/v = 0.064(3) 
is different from the exponents /3i and {32 ■ Within error 
bounds, Ijv = 0.93(1) is equal to l/vi and l/v 2 . All 
results of BP, SP, BPPR and SPPR are summarized in 
TableHl We can see that the critical exponents of BPPR 
are different from that of SPPR. The universality classes 
of BPPR and SPPR are different. 
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ln(L) 


3 4 5 6 7 8 9 
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Figure 9: Results of site percolation under the product rule: 
(A)log-log plot of r^(oo) — fc{L) with rc(oo) = 0.7634 and 
the exponent Xjvx — 0.42(2); (B)log-log plot of A(L), which 
increases at first and then decreases slowly; (C)log-log plot of 
Xr with the exponent 1 / 1/2 = 0.42(3); (D)log-log plot of xa 
with the exponent /32 = 0.008(4). 


Table I: Summary of critical points and critical exponents of 
BP, SP, BPPR and SPPR. rc, Pjv and 1/v of BPPR are taken 
from Ref.[l^. 


BP 

SP 

BPPR 

SPPR 

rc 0.5000(4) 

0.5927(3) 

0.5266(1) 

0.7634(2) 

l/u 0.75 

0.75(1) 

0.93(1) 

0.42(1) 

l/vi 0.75(1) 

0.75(1) 

0.99(5) 

0.42(2) 

l/i/2 0.74(1) 

0.74(1) 

0.94(3) 

0.42(3) 

00 

0 

0 

0.108(4) 

0.064(3) 

0.010(1) 

/3i 0.105(2) 

0.104(1) 

0.044(5) 


Pi 0.104(1) 

0.104(1) 

0.041(4) 

0.008(4) 


V. CONCLUSION 

We have studied the critical behaviors of the sizes and 
the size jumps of giant clusters on two-dimensional square 
lattice for both site and bond percolation under ran¬ 
dom and product rules respectively. Under the random 
rule, sites or bonds are added randomly into the lattice. 
We add the site or bond with smaller size product of 
two connected clusters into the lattice from two candi¬ 
dates picked randomly when the product rule is taken. 
The finite-size scaling forms of giant cluster size and size 
jumps are proposed and used to determine the transition 
points and critical exponents of the percolation transi¬ 
tions from Monte Carlo data. 

Our results show that the critical exponents of the size 
and size jump of random site percolation are equal to 
that of random bond percolation within error bounds of 
Monte Carlo data. As expected, the universality class 
of percolation is independent of site and bond under the 
random rule. The critical exponents of the size and size 
jump in the site percolation under the product rule are 
different from that in the random percolation and also 
the bond percolation under the product rule. The site 
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